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We study superfluid and Mott insulator phases of cold spin-1 Bose atoms with antiferromagnetic 
interactions in an optical lattice, including a usual polar condensate phase, a condensate of singlet 
pairs, a crystal spin nematic phase, and a spin singlet crystal phase. We suggest a possibility of 
exotic fractionalized phases of spinor BEC and discuss them in the language of topological defect 
condensation and Zi lattice gauge theory. 



Recent experiments on the condensation of Bose gases 
with internal degrees of freedom initiated a considerable 
amount of work on understanding the nature of spinor 
BEC. In particular Bose condensation in alkali atoms 
with nuclear spin 1 — 3/2 that have three low energy hy- 
pcrfine states and therefore behave as F = 1 bosons, have 
been a subject of active experimental Jl|-|| and theoret- 
ical research ]5[^To[. It was pointed out that the ground 
state of atoms with antiferromagnetic interactions in an 
optical trap is not a usual condensate of a macroscopic 
number of atoms in a single quantum state, but a com- 
plicated many body state of atoms arranged into a total 
singlet [R-ffo|- In this article we examine the problem 
of Bose F = 1 atoms with antiferromagnetic interac- 
tion (spin symmetric interactions are assumed through- 
out the paper) in optical lattices - arrays of microscopic 
potentials created by interfering laser beams 1 1 - lj]] . The 
dynamics of spinless bosonic atoms in arrays of optical 
wells may be described by a Bose-Hubbard model with 
the possibility of quantum phase transitions between in- 
sulating and superfluid phases induced by varying the 
properties of the laser light J2fj. Here, we propose sev- 
eral new phases that will appear for F = 1 bosons due to 
an interplay between spin and charge degrees of freedom: 
singlet pair condensate that only breaks charge symme- 
try, a spin nematic crystal phase that only breaks spin 
symmetry, and a novel "strong coupling pairing phase" 
that breaks both spin and charge symmetry but is dis- 
tinct from a simple polar BEC. These phases may be 
considered as "fragmented condensates" and have 

fractionalized topological excitations: half-vortices and 
7r-disclinations. Josephson type experiments may be used 
to distinguish among insulating and various superfluid 
phases. We also conjecture the possibility of fraction- 
alized phases in spinor Bose gases. Fractionalization is 
characterized by a topological order and may coexist with 
any broken continuous symmetry p2| . 

In what follows we consider the case of antiferromag- 
netic interaction between the atoms do > 02, where af 
is an s-wave scattering length in the F channel. In the 



individual wells we take wavefunctions that are superpo- 
sitions of the mean-field solutions of the interacting prob- 
lem in the absence of tunneling between the wells [p|jic|l 
\'*l>}N=f a '>l>(n)\N,n), with \N,n) = ^(n x al + %<4 + 
n z a\) N \Q). Here n is a real unit vector, are boson 
creation operators in the lowest vibrational state, |0) is 
the vacuum state, and M is a normalization factor. Two 
properties of these wavefunctions are important: (i) they 
satisfy a symmetry condition \N, — n) = (— ) N \N, n) [ ^o| ; 
(ii) in the limit of large N states that correspond to dif- 
ferent n's are orthogonal to each other. Hamiltonian of 
an array of identical optical wells with spin symmetric 
tunneling is given by |p3| , p4| 

H = 22 Hi + 2J ftij 
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Hij = -2tnnij {b\b 3 +b\bi) 



(1) 



where we defined the charge creation and annihilation op- 
erators that change the number of particles Ni, but not 
the direction of n,: b\\Ni,ni) = (Ni + 1) 1/2 |A^ + l,n 4 ); 



the number of particles in each well N = b\bi\ and 

d 

dm 
Par- 



the angular momentum operators = -in; x 
that describe a collective spin in well i 
ity condition on the wavefunctions implies constraint 
Ni + Li =even for all i, and for simplicity in (|l|) we 
chose to work in the grand canonical ensemble using 
a chemical potential /i. For a single well Hamilto- 
nian (Q) correctly gives a spin singlet ground state if 
the number of particles is even When opti- 

cal lattices are produced by lasers with wavelength A, 
they create an optical potential V(x) = Vq sin 2 (A:a;i). 
We can estimate parameters of the Hamiltonian ([!]) as 

g = ^E R ^^(^f/\ U = 2|^ (°0+2a 2 ) ( Va_ )3 /4 ; 

t = (AVoEn^Expi^i^y/^l - 2(^ 2 )}, where 
Er = H 2 k 2 /2M is the recoil energy. The average number 
of particles in each well will be related to the atomic gas 
density n as N = n (^) 3 (-ff-) 3 / 4 - The obvious parame- 
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ters to tune in experiments is N (e.g. by varying boson 
density) and the strength of periodic potential Vq (by 
tuning the laser intensity). When the latter is varied it 
affects tunneling exponentially and the interaction coeffi- 
cients only as a power law. In what follows, we therefore 
take U and g as constant and consider a phase diagram 
in the coordinates t and p.. 

There are two limits when construction of the phase 
diagram of (|l|) is simple. When g = there is no en- 
ergetic penalty for having high angular momentum (non 
zero spin) states of quantum rotors in (Q) and spin sym- 
metry is broken. For integer filling factors we have spin 
nematic insulator phases (NI) for small t and a super- 
fluid polar condensate phase (PC) for larger values of t. 
For fractional filling factors the system is always super- 
fluid. The other simple case is when g is large and states 
with odd number of atoms in individual wells are ener- 
getically costly: they are not allowed to have L = state 
of a rotor due to parity constraint and have a higher 
rotational energy. It is energetically favorable to have 
even numbers of atoms in each well arranged into sin- 
glet combinations, which can be interpreted as binding 
of atoms into singlet pairs with a pair binding energy 
E P (Q) = 2E(Q + 1) - E(Q + 2) - E(Q) = g. In the 
limit when g — > oo the crystal phases are possible for 
even numbers of atoms per well only and correspond to 
spin singlet insulators (SSI). The superfluid phase is a 
condensate of singlet pairs (SSC), in which tunneling of 
individual atoms between the wells is suppressed and only 
singlet pairs are delocalized. The origin of pairing in this 
case is not the attraction between individual bosons, but 
a singlet formation on the scale of individual wells. This 
is reminiscent of the "attraction from repulsion" mech- 
anism of electron pairing proposed by Chakravarty and 
Kivelson for high Tc cuprates, Cqq, and polyacenes [gs) . 
In the case of finite g we expect Mott crystal phases for 
all integer filling factors. For even N's and small t the ro- 
tors kinetic energy dominates and we have singlet crystal 
phases (SSI). When t is increased the system goes to a 
spin nematic insulator (NI) phase that has admixture of 
L 7^ in the ground state. For odd iV's when all charge 
fluctuations are frozen there is always at least L = 1 on 
each site, so we expect the system to be analogous to 
spin 1 lattice model with a broken spin symmetry. The 
superfluid phase will be spin singlet fluid of pairs (SSC) 
for small t, and a polar condensate (PC) for larger t's. 

Simple Josephson type experiment may be used to dis- 
tinguish superfluid phases suggested above. We imagine 
spinor atoms in an optical lattice in the presense of grav- 
ity. Gravity acts as voltage in a conventional Josephson 
junction arrays and gives rise to Josephson oscillations 
jl9| . Josephson oscillations may be measured selectively 
in spin (1), or measuring all the particles regardless of 
their spin (2). In the case of single atom condensate 
(PC) we expect oscillations for both kinds of measure- 
ments with frequent hujj = mgX/2. A singlet pair con- 



densate will only give oscillations in channel (2) with fre- 
quency 2huj = mg\. Insulating phases will not show 
any Josephson oscillations, and nematic vs spin singlet 
may be distinguished further through different correla- 
tions for particles of different spin. 
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FIG. 1. Phase diagram of the Hamiltonian (|lj) for different 
values of parameter g. 

Order parameters for various phases suggested above 
and shown on Figure [l] may be conveniently discussed 
using a Euclidean space-time action. We implement 
a constraint N% + L, =even using a projection opera- 
tor Pi = |E Ji=±1 e'5 (1 ^ )(,Vi+li) @. To calculate 
the partition function Z = Tr[e~"' H P] we divide imag- 
inary time (3 into M slices of length e = j3/M and af- 
ter a few standard manipulations that include Hubbard- 
Stratanovich transformation of the last term in (|l]) we 
find 

S = — Jr r ,(7 rr >C0S(j) rr > — Jrr'COs(2(f) rr f) 
rr' rr' 

- ^2 Jr r >(7rr>n. r n r > - ^ J^Q° 6 Q"? (2) 
rr' rr' 

Here summation goes over sites r = (z,t) in the space- 
time lattice, <fi r is the phase variable conjugate to N r , 
<fi rr * = <fi r — (fi r >, Qr b — n^n b r — S ab /3 is a nematic order 
parameter, and a rr i = ±1 is an Ising field that lives on 
the links rather than sites. The coupling constants are 



j2c _ j2s 
r,r±r r,r±r 



J. 



r.r±f ~ ( e d) 
2c, 2s 



0, J 



r,r±{x,y,z} 



Jr,r±{x,y,z} ~ e ^lxl' 

— eJ/4, with x being 
the saddle point value of the Hubbard-Stratanovich field. 
The usual periodic boundary conditions for 0, n, and a 
are assumed at r = and r = (3. In writing (0) we 
omitted the Berry phase terms that are important for 
quantitative calculations of the phase diagram but not 
for the symmetry arguments discussed below. 

Action (^) has a Z2 gauge symmetry <fi r — » cfi r + ^(1 — 



6 r ), n r 



r,r+a t r-fa? 



where e r = ±1. 



Such a symmetry has been pointed out earlier in [ p_0| , 
where it was observed that the physical order parameter 
is a complex vector d = e'^n (not <fi and n separately) 



that has a symmetry of e 1 ' 



J<t> 



Equa- 



tion (Eh is the simplest action consistent with the charge 
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U(l), spin 50(3), and gauge Z 2 symmetries of the model. 
Another term allowed by the Z 2 symmetry is the ana- 
logue of the Maxwell terms for the lattice gauge models 
S a = —K J2nTln a ijj where summation goes over pla- 
quettes in d + 1 dimensional lattice. This term may be 
generated by integrating out the high energy degrees of 
freedom. 

The existence of the local symmetry imposes impor- 
tant constraints on the possible order parameters and 
symmetry breaking states of the system. Only order pa- 
rameters that are gauge invariant may acquire expecta- 
tion values. For example, the expectation values (n) ^ 
or (e 1 ^) ^ are not allowed in the models described 
by (||), since n, and e l ^ i are not invariant under Z 2 
gauge transformations. Physically this restriction arises 
from the fact that the wavefunction ^(n) should be an 
even or odd function of n for N even or odd respec- 
tively. Below we consider several examples of the order 
parameters that are gauge invariant and review broken 
symmetry phases that they lead to (see Figure |l|). Po- 
lar BEC phase (PC) has (d) = (e^ n) ^ 0. It breaks 
both charge and spin symmetries and is the phase where 
atoms are condensed directly. Nematic Insulator phase 
(NI) has {Qab) — (n a nb — S a b/3) ^ and breaks only 
the spin SO(3) symmetry. This is a phase where not 
the atoms but certain particle-hole composites are con- 
densed. Spin singlet condensate (SSC) (e 12 ^) ^ breaks 
the U(l) symmetry but not 50(3). It corresponds to 
a condensate of singlet pairs of atoms. Strong coupling 
pairing phase (SCP) is a phase where both (e l2< ^) ^ 
and (Q ab ) ^ 0. This phase breaks both U{\) and 50(3) 
symmetries, however, it is fundamentally different from 
the polar condensate phase. Transition between phases 
PC and SCP should be Ising type transition. Similar 
phase has been recently suggested in the context of bi- 
layer quantum Hall systems and called "strong coupling 
pairing phase" pq| . Finally, spin singlet insulator (SSI) 
is a phase where no field has expectation value, so neither 
U(l) nor SO (3) symmetries are broken. We note that for 
a finite number of traps, when spontaneous symmetry 
breaking is not possible, one can observe signatures of 
the transitions by measuring fluctuations in the number 
of atoms, atom pairs, and spins JHJ. 

It is interesting to point out that phases NI, SSC, and 
SCP correspond to the fragmented BEC discussed by P. 
Nozieres and D. Saint James. The part emphasized in 
plf is that fragmented state has no macroscopic popula- 
tion of k=0 state for the bosons. Another interesting as- 
pect of these states is that they achieve fractionalization 
of topological objects. The unfragmented PC phase has 
topological excitations that are composites of 1/2 vortex 
in the charge sector and 7r-disclination in the spin sector. 
Here 1/2 charge vortex is a topological defect around 
which the phase 4> winds by ir and e 1 ^ changes sign, ir- 
disclination in the spin sector is introduced as a topo- 
logical defect around which n changes sign (vortex like 



objects in n are called merons, so 7r-disclination is 1/2 of a 
meron). Both 1/2 vortex and 7r-disclination are pointlike 
objects in two dimensions and lines in three dimensions. 
No individual 1/2 charge vortices or 7r-disclinations are 
allowed in the PC phase. The spin and charge parts of 
the vortex are always glued together. A condensate of 
singlet pairs only, the SSC phase, has 1/2 charge vor- 
tices as separate excitations, but no spin vortices; and 
the nematic phase NI has ir disclinations and no charge 
vortices. A condensate of S=0 pairs and S=2 pairs, a 
SCP phase, will have 1/2 charge vortices and ir disclina- 
tions as separate excitations. 

Establishing the nature of the broken symmetry does 
not fully characterize models with gauge symmetry. An- 
other important aspect of the system described by (||) 
is the possibility of confining and deconfining phases of 
the Z 2 gauge theory in dimension d + 1 > 3 |2^,|27],^8| . 
Difference between the confining and deconfining phases 
relies on the concept of topological order (29) and has im- 
portant implications on the nature of excitations in the 
system. Topological order may coexist with any true long 
range order, so any of the phases reviewed earlier may be 
confining or deconfining. For the pure Z 2 gauge models 
the confining and the deconfining phases may be distin- 
guished using Wilson loops, which are characterized by 
an area law for the confining phase and by a perimeter 
law for the deconfining phase. With matter field present 
Wilson loops may no longer be used to discriminate be- 
tween the two phases ]2^], however distinction between 
them survives, and there is a phase transition between 
the two. A simple way to understand the transition is to 
think of it as condensation of visons, topological excita- 
tions of the Z 2 gauge theory that describe frustrated pla- 
quettes Yin a rr' = — 1- Particles that carry a Z 2 charge 
are frustrated when traveling around a vison, so when 
visons are condensed such particles may not propagate 
coherently and we have a confining phase. When visons 
are not condensed, they are finite energy topological ex- 
citations and we have a deconfined phase. Presense of 
finite energy visons leads to an additional degeneracy of 
the gauge models in a deconfined state on manifolds with 
non-trivial degeneracy. Systems with Z 2 gauge symme- 
tries in condensed matter have already been discussed 
by Lammert et. at [ [30[ in connection with liquid crystals 
and Senthil and Fisher in connection with high tempera- 
ture superconductors |22]| . For our system the discussion 
above implies that all the broken symmetry states , as 
well as a state SSI, where neither U(l) nor SO (3) sym- 
metries are broken, come in two varieties: confining and 
deconfining. Detailed comparison between confining and 
deconfining versions of various will be given in |23j . Here 
we only note that some of the most striking implications 
of deconfinement appear for the SSI phase. The deconfin- 
ing SSI* phase allows charge b = e 1 ^ and spin n carrying 
excitations propagate independently pl| |. 

Another perspective on the nature of deconfining 
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phases comes from considering phase transitions as con- 
densation of topological defects. Discussion will be given 
in the case of two spatial dimensions, where condensation 
of point-like topological objects (vortices, 7r-disclinations, 
merons, and etc.) drive transitions into the phases of re- 
stored symmetry. We expect, however, that most of the 
conclusions reached for d = 2 will also hold in d — 3. 
Let us consider a transition between the PC phase that 
breaks charge and spin symmetry and NI phase that 
breaks only the spin symmetry. In the phase PC we have 
topological objects that are 1/2 vortex 7r-disclination 
composites. Half-vortex and half-meron parts may come 
with different signs, so we have 4 types of these com- 
posites (±1/2, ±1/2), where the first and the second 
numbers stand for the charge and the spin parts respec- 
tively. To go to a phase that has broken spin sym- 
metry and restored charge symmetry one expects that 
we need to condense topological objects in the charge 
sector, with topological objects in the spin sector re- 
maining at finite energy. This can be achieved by con- 
densing composites (1/2, 1/2) + (1/2, -1/2) = (1,0) and 
(-1/2,1/2) + (-1/2,-1/2) = (-1,0), but having indi- 
vidual topological defects (±1/2, ±1/2) gapped. Such 
transition would correspond to condensation of integer 
vortices and would be equivalent to order-disorder tran- 
sition for the charge order parameter e 1 ^ in the usual 
Bose-Hubbard model. An important observation is that 
the insulating phase obtained this way is a fractional- 
ized phase NI*! It does not have a condensate of mini- 
mal topological object in the charge sector (1/2 vortex), 
but only a condensate of integer vortices. So, only the 
PC — NI* transition may be a continuous second order 
transition, but PC — NI must be a first order transition. 
The existence of fractionalized phases may, therefore, be 
established through study of phase transitions pb|| . We 
note, however, that starting from a strong coupling pair- 
ing phase SCP (that has unbound 1/2- vortices and ir- 
disclinations) we will have SCP — NI as a second order 
transition. 

To summarize, we reviewed several novel phenomena 
that will take place for F — 1 Bose gases with antifcrro- 
magnetic interactions. We showed that in optical lattices 
they will have several distinct superfluid and Mott insu- 
lator phases depending on the density and parameters of 
the system. Another possibility suggested in this paper 
is a class of fractionalized phases of spinor BEC, char- 
acterized by the topological order that may coexist with 
any true long range order in the system. Quantitative 
analysis of experimental systems will be given in subse- 
quent publications ]23| ] . Discussion in this paper may be 
generalized to systems with higher hyperfine spins. We 
also expect some of the results to be applicable to triplet 
superconductors. 
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